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Hadron Masses and Quark Condensate from Overlap Fermions 
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We present results on hadron masses and quark condensate from Neuberger's overlap fermion. The scaling and 
chiral properties and finite volume effects from this new Dirac operator are studied. We find that the generalized 
Gell-Mann-Oakes-Renner relation is well satisfied down to the physical u and d quark mass range. We find that 
in the range of the lattice spacing we consider, the n and p masses at a fixed m^-/m p ratio have weak 0(a 2 ) 
dependence. 



The recent advance in chiral fermion formula- 
tion which satisfies Ginsparg- Wilson relation has 
a great promise in implementing chiral fermion 
for lattice QCD at finite lattice spacing. It is 
shown to have exact chiral symmetry on the lat- 
tice [Q and it has no order a artifacts Q. Neu- 
berger's Dirac operator |l| derived from the over- 
lap formalism has a compact form in four dimen- 
sions which involves a matrix sign function 



with Q =/= 0, we find the critical slowing down 
is much milder than that of the Wilson fermion 
and there are no exceptional configurations. The 
critical slowing down sets in quite abruptly after 
[ia = 0.003 which is already at the physical u and 
d masses. 

It is shown |J that the generalized Gell-Mann- 
Oakes-Renner (GOR) relation 



D=~[l+ t i+(l- t x) l5 e(H)]. 



(1) 



d 4 a;(7r(a;)7r(0)) = 2(#*) 



(2) 



In this talk, I present some preliminary re- 
sults from our numerical implementation of the 
Neuberger fermion. We adopt the optimal ratio- 
nal approximation of the matrix sign function [Q 
with 12 terms in the polynomials. The small- 
est 10 to 20 eigenvalues of H 2 are projected out 
for exact evaluation of the sign function for these 
eigenstates. We use multi-mass conjugate gradi- 
ent as the matrix solver for both the inner and 
outer loops. With residuals at 10 -7 , the inner 
loop takes ~ 200 iterations and the outer loop 
takes ~ 100 iterations. We check the unitarity 
of the matrix V — jse(H). For Vx — b, we find 
\x*x — b'b\ ~ 10~ 9 . Even for topological sectors 
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with tt(x) being the pion interpolation field, is 
satisfied for each quark mass and volume, config- 
uration by configuration. We utilize this relation 
as a check of our numerical implementation of the 
Neuberger operator. We find that for the lattices 
we considered (6 3 X 12, p = 5.7, 5.85; 8 3 X 16, P = 
5.85, and 10 3 x 20, P = 5.85,6.0) the GOR rela- 
tion is satisfied very well (to within 1%) all the 
way down to the smallest mass fia = 0.0001 for 
the Q — sector. For the Q / sector, the 
presence of zero modes demands higher precision 
for the approximation of e(H). For example, we 
show in Fig. [I] the ratio of the right to left side 
of Eq. (H) for a configuration with topology on 
the 6 3 x 12 lattice at p = 5.7 as a function of the 
quark mass. When only 10 smallest eigenmodes 
are projected, we see that the ratio deviates from 
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one for small quark masses. The situation is con- 
siderably improved when 20 smallest eigenmodes 
are included. The situation is better than the 
domain-wall fermion case when the size of the 
fifth dimension is limited to L s = 10 to 48 HI. 



0.95 
0.9 



20 eigenmodes 
10 eigenmodes 



Figure 1. Ratio of the right to left side of Eq. 
(^) for a configuration with topology. <>/• indi- 
cates the case with projection of 10/20 smallest 
eigenmodes. 

We also calculate the quark condensate (W^) 
with 3 to 6 Zi noises for each configuration. For 
small quark mass, it has the form 



( ^, ) = M> +Co + Cl/i . 



(3) 



The singular term which is due to the zero modes 
in the configurations with topology (Q ^ 0) is 
specific to the quenched approximation. It will 
be suppressed when the determinant is included 
in the dynamical fermion case. We see this clearly 
in the following figure which is first seen with the 
domain- wall fermion ||. A fit to the formula in 
Eq. (||) is given in Fig. ^. We see that Co is non- 
zero. The standard definition of the quark chiral 
condensate entails the extrapolation of Cq to the 
infinite volume before taking the masslcss limit. 
Another way is to consider the finite-size scaling 
@. When the size of the lattice is much smaller 
than the pion Compton wavelength, i. e. L <C 
1/rrin, the (^^>) is proportional to [iY?V for small 
masses besides the - term due to quenching. 
From this, the infinite volume condensate £ can 
be extracted. We plot in Fig. || (^^)a 3 / \xa vs fia 
in the Q — sector for 3 lattice volumes (6 3 x 
12, 8 3 x 16, and 10 3 x 20) at [3 = 5.85. We see 
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Figure 2. Quark condensate as a function of the 
quark mass. This is done on 50 6 3 x 12 lattices 
at j3 = 5.7. 
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Figure 3. (^"3>)a 3 / fia vs fia in the Q = sector 
for 3 lattice volumes at j3 = 5.85. 



that they are quite flat which indicates that the 
condensate is indeed proportional to fi and we 
also see that they increase with volume. 

We have calculated the ir,p and nucleon 
masses. A typical result on the 8 3 x 16 lattice 
at (3 = 5.85 is given in Fig. |J. We see the fi- 
nite volume effect on the nucleon mass when fia 
is smaller than ~ 0.15. To see the behavior of 
pion masses near the chiral limit, we plot mja 2 
as a function of fia in Fig. || for three lattices with 
about the same physical volume. It appears that 
there might be a ^/Jla behavior in the very small 
\xa region which we will explore further. When 
we project only 10 smallest eigenmodes in the ap- 
proximation for the sign function in the 6 3 x 12 
case, we see that m^a 2 tends to a finite value as 
fia — > 0. This implies a residual mass due to the 
poor approximation of e(H), a behavior similar 
to that observed in the domain- wall fermion with 
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10 J x 20 with )3 = 6.00 
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Figure 4. Masses of n,p and N on the 8 3 x 16 
lattice at j3 = 5.85 are plotted vs \xa. 



finite L s . Finally, we check scaling. We plot in 
Fig. H Tr/y/a and p/y/a vs aa 2 where a is the 
string tension from which the lattice spacings are 
determined Q . It is known that the overlap oper- 
ator does not have 0(a) artifacts. Now it appears 
that the 0(a 2 ) errors are small. 
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Figure 5. m^a 2 are plotted as a function of \ia 
for three lattices with about the same physical 
volume. 
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Figure 6. n / '■y/a and p/y/cr are plotted vs aa 2 for 
three lattices with similar physical volumes. 



